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Abstract 

In this paper, we provide a simple, "generic" interpretation of multifractal scaling laws and 
multiplicative cascade process paradigms in terms of volatility correlations. We show that in 
this context 1/f power spectra, as recently observed in Ref. p0[ , naturally emerge. We then 
propose a simple solvable "stochastic volatility" model for return fluctuations. This model is able 
to reproduce most of recent empirical findings concerning financial time series: no correlation 
between price variations, long-range volatility correlations and multifractal statistics. Moreover, 
its extension to a multivariate context, in order to model portfolio behavior, is very natural. 
Comparisons to real data and other models proposed elsewhere are provided. 



1 Introduction 

As shown by most recent empirical studies on huge amount of data, the market price changes 
are characterized by several "universal" features Jj], |^]: price increments are not correlated, 
volatilities are strongly (power-law) correlated and price increment probability density function 
(pdf) shapes depend on the time scale. From quasi Gaussian at rather large time scales, these 
pdf are characterized by fat tails at fine scales. Many authors in the recently emerged field 
of "econophysics" Jj], |^, |j| as well as in classical empirical finance, aim at proposing simple, 
discrete or continuous time models that are able to account for these observations. Among all 
the proposed models, one can distinguish several streams, from the simplest Brownian process, 
that constitutes the main tool used by practitians, to the class of "heteroskedastic" nonlinear 
processes as proposed in Refs. f| [?|. To account for the letpokurtic nature of the small scale 
pdf, Mandelbrot Q and Fama [|| proposed the Levy stable paradigm that has been recently 
improved in the "truncated Levy" version [|J §]• More recently, an interesting comparison 
between market price variations and the fluctuations of the fluid velocity field in fully developed 
turbulence has been suggested || . Besides the real pertinence of such an analogy that has been 
widely commented JTT| , |l2| , , this work opens the door to another important paradigm to model 
financial time series, namely multifractal processes. The multifractal processes^, and the deeply 
connected mathematics of large deviations and multiplicative cascades, are well known to be 
useful to describe the intermittent nature of fully developed turbulence [ pl| . Recent empirical 

people sometimes refer to "multi-amne" processes or processes that display "multi-scaling" 
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findings jH], [L3| [15], [16| suggest that in finance, this framework is also likely to be pertinent as 
far as the time scale dependence of the statistical properties of price variations is concerned. 

Our purpose in this paper is twofold. First, we make a brief review of multifractals in order 
to specify what is a multifractal process. We try to provide several complementary points of view 
and to understand what are the main ingredients for "multi-scaling" . We also comment about 
the criticisms raised by several authors about multifractality in finance. Our second goal is to 
propose a simple multifractal "stochastic volatility" model that captures very well all the above 
mcntionned features of financial fluctuations. This model, that has been originally introduced in 
Ref. Pq|, is compared to real data and some models proposed elsewhere. We discuss its possible 
multivariate extension in order to use it in management applications. The paper is organized 
as follows. The review on multifractal processes is made in section 2. We introduce notations, 
the related notions of multi-scaling, scalc-invariance, cascade process and self-similarity kernel. 
We illustrate our purpose using empirical estimates for some high frequency financial data. In 
section 3 we review some findings of Ref. | flo| concerning the magnitude correlations for cascade 
models and suggest a link with 1// processes as recently observed in financial time series. In 
section 4 we introduce the multifractal random walk defined in Ref. |i| as a stochastic volatility 
model. We discuss its main properties and propose a natural multivariate generalization. Our 
discussion is illustrated by numerical simulations. In section 5 we propose estimators for the 
few parameters of our model and compute them for some intraday and daily time series. In 
section 6 we discuss some related works about multifractality in finance. Conclusions and some 
prospects are reported in section 7. 



2 Multifractal processes and cascade models 

In this section we briefly discuss the related notions of multifractality and multiplicative cas- 
cade. Most of the ideas and concepts that we recall below have been introduced in the field 
of fully developed turbulence where people aim at accounting for the so-called "intermittency 
phenomenon" (for a review of this subject see e.g., (H]). 



2.1 Multifractality of financial time series 

Let us consider the variations of a stochastic process X(t) at a time scale I. For that purpose, 
one can consider the increments of the process, §iX(t) = X(t + I) — X(t) or more generally its 
wavelet transform |2^, [2^, Q 

T{t,l) = l~ l j ' i>(t^£\x{t')dt' 

where tp(t) is the so-called analyzing wavelet, i.e, a function well localized in both Fourier and 
direct spaces^. Let us denote M(q,l) the order q absolute moment of SiB(t) or T(t,l), (in this 
paper E(.) will be used for the mathematical expectation and we will always suppose that the 
considered processes has stationary increments) 

M(qJ) = E(\6 l X(t)n. (1) 

We will say that the process is scale-invariant, if the scale behavior of the absolute moment 
M(q, I) is a power law. Let us call ( q the exponent of this power law, i.e., 

M(q,t)~C q fr, (2) 



2 One nice property of wavelet transform is that it can be inverted, i.e., one can recover the original signal from its 
wavelet coefficients. Another interesting feature is that there exist orthonormal wavelet bases. Such bases are very 
useful for signal synthesis and modelling, as it is illustrated for cascade processes in Ref. |33] 
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Figure 1: Multifractal Analysis of the intraday future S&P500 index over the period 
1988-1999. (a) Plot of the original index time-series. The analyzed time-series is the detrended 
and de-seasonalized logarithm of this series, (b) Log-log plots of M(q, I) versus I for q = 1, 2, 3, 4, 5. 
The time scales I range from 10 minutes to 1 year, (c) log 2 (M(g, l)/M(l, l) q ) for q = 2, 3, 4, 5. Such 
plots should be horizontal for a process that is not multifractal. (d) C q spectrum for the S&P 500 
fluctuations. The plot in the inset is the parabolic nonlinear part of ( q . 
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where C q is a prefactor that will be interpreted below. The process is called monofractal if 
C q is a linear fonction of q and multifractal if Cq is nonlinear. Note that, from the concavity 
of the moments of a random variable, it is easy to show that Cq, as defined from the scaling 
behavior (j^) in the limit / — > + , is necessarily a convex fonction of q. The same argument leads 
to the conclusion that such scaling behavior with a nonlinear Cq cannot hold for all scales 
Thus, for a multifractal process there exists at least one characteristic time T (hereafter referred 
to as the integral time) above which the behavior (^|) is no longer valid. Multifractality has 
been introduced in the context of fully developed turbulence in order to describe the spatial 
fluctuations of the fluid velocity at very high Reynolds number pll. As suggested by recent 



studies [|9| |10|, |13j |15|, |16fl , multifractality is likely to be a pertinent concept to account for the 
prices fluctuations in financial time-series. This is illustrated in Fig. [I] where the Cq function 
is estimated for the future S&P500 index over the period 1988-1999. The original intraday 
time-series has been sampled at a 10 mn rate (Fig. |l|(a)) in order to obtain equi-sampled data. 
We consider the associated continuously compounded return time-series, (i.e., the logarithm of 
the index value) that has been detrended and de-seasonalizedf]. The Cq spectrum in Fig. |l](c) 
is obtained using linear regression fit of "log-log" representations of the behavior of the g-th 
order moment versus the time scale as illustrated in Fig. |l](b). In this figure, the scales span an 
interval from 10 minutes to approximately 1 year. Moment estimates at larger time scales are 
very poor because of the finite size of the overall record. From the linear behavior of such curves, 
one clearly sees that the scale- invariance hypothesis is satisfied over around 3 decades. In Fig. 
|](b) we have plotted log 2 m{i,i)i versus l°g2(0- The fact that such plots are not constant reflects 
the nonlinearity the Cq spectrum. The future S&P500 can thus be considered, at least at this 
description level, as a multifractal signal. Let us notice that we have computed, in Fig. |l|(d), the 
Cq values for g-th order moments that include negative values of q. This can be achieved using 
a wavelet based technique that has been introduced in Refs. [^5[ p6[ [27], |28|. This spectrum 
turns out to be well fitted by a parabolic shape Cq — 0.53g — 0.015g . The non linear parabolic 
component of Cq has been plotted in the inset of Fig. [l](d). 



2.2 Multifractal processes, self-similar processes and multiplicative 
cascades. 

Multifractality (in the sense defined above) is a notion that is often related to an underlying 
multiplicative cascading process. In the context of deterministic functions the situation is rather 
clear since the analyticity of the Cq spectrum is deeply connected to the self-similarity properties 
of the function [ p6| p7| , |28| , p9[ . Roughly speaking, a function is self-similar if it can be written as 
a multiplicative cascade in an appropriate space-scale (or time-scale) representation [^7], |2^| . 
In that context, the so-called multifractal formalism is valid, i.e., one can relate the Cq spectrum 
to the D(h) singularity spectrum that provides information about the statistical distribution 
of singularity (Holder) exponents h. The things are somehow more complex for stochastic 
processes. One of the goals of this paper is to provide some simple elements about this subject. 

In the mathematical literature, a process X(t) is called self-similar of exponent H if VA > 0, 
\~ H X(Xt) is the same process as X{t). According to this definition, the Brownian motion is 
self-similar with an exponent H = 1/2. This definition is however too restrictive for our purpose 
since it excludes multifractal processes. Indeed, let us consider Pi(5X) the probability density 
function (pdf) of (^X(i)Q. If X(t) is self-similar with an exponent H, then it is easy to prove 

3 The amplitude of the return variations in each intraday period is normalized according to the estimated U-shaped 
intraday r.m.s. 

4 Note that from stationarity of the increments, the law of 5iX(t) is the same as the law of X(l) if one assumes 
that X(0) = 0. 
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that 

Pi{5X) = X H P M {X H 8X) . (3) 
Then, the moments at scale I and L — XI are related by 

L 



M(q,l) = C q ( T ) , (4) 



with C q = M(q, L). Thus one has a "monofractal" process with £ g = qH. In order to account for 
multifractality, one has to generalize this classical definition of self-similarity. This can be done 
by introducing a weaker notion, as originally proposed in the field of fully developed turbulence 
by B. Castaing and co-authors |3lj]. According to Castaing's definition of self-similarity, a 
process is self-similar if the increment pdf's at scales I and L = XI (A > 1) are related by the 
relationship j3l], |32| : 

Pi {8X) = J G hL {u)e- u P L (e- u 8X)du , (5) 

where the self- similarity kernel G^l depends only on l/L. Let us note that this definition 
generalizes Eq. (0) that corresponds to the "trivial" case G/ i l(u) = 8(u — H\n(l/L)). This 
equation basically states that the pdf Pi can be obtained through a "geometrical convolution" 
between the kernel Gi t L and the pdf Pi,. A simple argument shows that the logarithm of the 
Fourier transform of the kernel Gi l can be written as F hL {k) = \xiGi, L {k) = F(k)\n(l/L) Q. 
Thus, from Eq. (@), one can easily show that the q order absolute moments at scales I and L 
are related by: 

l \ F(-iq) 



M(g,l) = G l<L (-iq)M(q,L) = M(q, L) I -J , (6) 

and then C q = M(q, L) and C, q — F(~iq). A nonlinear (q spectrum implies that F is nonlinear 
and thus that G is different from a Dirac delta function^. For example, the simplest non linear 
case is the so-called log-normal model that corresponds to a parabolic ( q function and thus to 
a function G that is Gaussian. 

The equation (||) can be interpreted as follows: the pdf at scale I, Pi is written as a weighted 
superposition of the rescaled versions of the pdf at scale L, Pl, the self-similarity kernel Gi } l 
being the associated distribution of weights. In the case of a monofractal process as described 
by Eq. (||), a single value of u is sufficient in the equation (|J) since Pi and Pl have the same 
shape and differ only by the scale factor e _u = (l/L) H = X H . This explains the Dirac function 
for the kernel G. This situation can be easily generalized by considering other shapes for the 
kernel G/.l. In that case, the shapes of the pdf P; across scales are no longer the same: when 
going to small scales, fat tails emerge and the pdf become strongly leptokurtic (see fiefs. ||, |3l| 
or Fig. ||). 

Let us now make the link with multiplicative cascades. This can be easily done if one consider 
discrete scales /„ = 2~ n L. Let us suppose that the local variation of the process 8i n X at scale 
/„ is obtained from the variation at scale L as 

5 ln X(t) = (f[W^J 5 L X(t) (7) 

where Wi are i.i.d. random positive factors. This is the cascade paradigm. Realizations of such 
processes can be constructed using orthonormal wavelet bases as discussed in Ref. |33] . If one 



5 It essentially results from the fact that Gi t L depends only on l/L and satisfies the semi-group composition law 
Gi l ,i 3 = Gi lt i 2 * G; 2 ,; 3 where h < h < ^3 and * is the convolution product fnl, [sif . 

6 Note that from the above mentionned semi-group property, the Levy theorem J30| implies that G is necessarily 
an infinitely divisible law 
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defines the magnitude u)(t, I) at time t and scale I as the logarithm of "local volatility" |iofl : 

c(U) = iln(|<W)| 2 ), (8) 

then the previous cascade equation becomes a simple random walk equation, at fixed time t, 
versus the logarithm of scales: 

u(t, l n+ i) = u(t, l n ) + ln(W„+i) . 

If the noise InWi is normal N(fi,X 2 ), the pdf of u>, Pi(ui), thus satisfies a simple diffusion 
equation with a Gaussian kernel: 

Fi>)= (N(jjL,X 2 ) m *p L )(u) (9) 

where * is the convolution product. Going back to the original variable SX, the previous 
equation corresponds exactly to Castaing's formulation of self-similarity (|J) with the log-normal 
propagator: 

G lniL = N(^X 2 )* n = N(n^nX 2 ) . 

Conversely, let us consider a process that satisfies Castaing's equation with a normal kernel G. 
This means that one can write, 

SiX(t) = WS 2 iX(t) (10) 

where = means the equality in law of the two random variables and W is a log-normal random 
variable which mean /i and variance A 2 do not depend on I. By iterating this equation n times, 
one thus recover, at least heuristically, the cascade equation (||). Thus, the cascade picture across 
scales, constitutes a kind of paradigm of non-trivial self-similar processes. As explained in Ref. 
lILsf , the problem with such processes is that they involve representations (e.g., orthonormal 
wavelet bases) that are constructed on a discrete set of scales (e.g., dyadic scales l n = 2~ n ) and 
in turn cannot be invariant under continuous scale dilations. 



3 Magnitude correlations and 1/f spectra 

We have seen in the previous section that multifractality can be interpreted as a diffusion of 
the magnitude of the variations of the return from large time scales to small time scales. In the 
financial framework, magnitudes at all scales are nothing but a logarithmic representation of 
local volatilities. In this section we would like to address the problem of volatility correlations. 
The "heteroskedastic" nature of financial time-series is now a well established empirical fact. 
Volatility possesses long-range positive correlations: periods of strong activity alternate with 
quiet periods. A lot of models have been proposed to account for this phenomenon from the 
famous GARCH models to various stochastic volatility models. Let us proceed with the multi- 
fractal and cascade picture and study what kind of correlations are associated to these models. 
This problem has already been considered by Arneodo, Muzy and Sornette in Ref. (see 
also Refs |33], M). These authors have shown that a log- normal cascade model on the dyadic 
tree associated to the orthonormal wavelet representation leads naturally to magnitude corre- 
lation functions G u {1,t) — Cov(ui(t,T),uj(t + i,r)) that behave as — X 2 \n(l/T) for T > I > r. 
This behavior has been shown to provide good fits of the empirical estimates of the correlation 
functions from real data p0| . In Fig. || is reported the magnitude correlation function C u (l, r) 
(we choose r = 10 min) of the S&P500 time series studied in Fig. [l| One can see that, when 
plotted versus the logarithm of the time lag, \n(l), the correlation function decreases linearly 
with a slope A 2 ~ 0.025. The intercept of such straight line provides an estimator of the integral 
time T that is, in our case, approximately T ~ 3 years (note that because we get an estimate 
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Figure 2: Magnitude correlation function of the S&P 500 future, (a) C W (Z, r) versus Z for 
r = 10 min. The solid line represents a fit according to the cascade model logarithmic expression, 
(b) Cu{1,t) versus ln(Z). The cascade model predicts a linear behavior that crosses the y-axis at 
ln(Z) = ln(T). The small scale cross-over is due to the smoothing window used to estimate the local 
magnitude uj(t,r). 
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of ln(T) the error on the value of T is very large). We have checked that those results are 
stable when changing the reference time r for return calculation. As it will be illustrated in 
section 5, for various financial time series, the "cascade ansatz" is very pertinent to describe the 
volatility correlations. Let us notice that the very slow (logarithmic) decrease of the correlation 
functions for time lags below the integral time T, is reminiscent of the ultrametric nature of 
the tree naturally associated to the time-scale (or time-frequency) representation |n| [33| |34| . 
Moreover, let us remark that, as far as power spectrum is concerned, Gaussian processes with 
such correlation functions can be seen as "1//" processes. Indeed, if the correlation function is 
given by the above expression, the power spectrum can be shown to reduce to 

S(f) = 2A 2 /" 1 / f x- 1 sm(x)dx . (11) 
Jo 

In the high frequency limit / — > +oo, we then have S(f) ~ A 2 7r/ _1 . Another intuitive way 
to understand this property, comes from the fact that the logarithmic decay of the correlation 
function can be understood as the limit H — > in the power-law correlation function k~ 2H of a 
fractional gaussian noise of exponent H. This property will be explicitely used in the discussion 
of section 6.3. Let us finally remark that 1/f spectra have been observed in a wide range 
of applications Recently, Bonanno et al. Q suggested the possible pertinence of such 

processes to account for the fluctuations of the number of trades of different stocks. 



4 A simple solvable multifractal model 

As emphasized previously, multiplicative cascade models represent the paradigm of multifrac- 
tal processes in that they contain the main ingredient leading to multifractality, i.e, the scale 
evolution of the magnitudes, from coarse to fine scales, is a random walk. Besides the problems 
of continuous scale invariance and stationarity of standard hierarchical constructions of such 
processes, they cannot be formulated using a stochastic evolution equation as one would expect 
for a model for financial time series. In this section we propose a "stochastic volatility" model 
that has been introduced in Ref. |l8| , that does not possess any of these drawbacks: it has 
stationary increments, it has log-normal multifractal properties and is invariant under contin- 
uous dilations. The key idea underlying this model is that the stochastic volatility possesses, 
as for cascading processes, a "1//" spectrum, or, more precisely, a correlation function with a 
logarithmic behavior. 



4.1 The multifractal random walk 

Let us briefly recall the construction of the multifractal random walk (MRW) proposed in |is| . 
A discretized version of the model X At (using a time discretization step At) is built by adding 
up tj At random variables : 

t/At 

X A t(t)=J2 SX ^[k}, 

k=l 

where the process {c>A'At[fc]}/c is a noise whose variance is stochastic, i.e., 

6X At [k] = e At [k]e u ^ , (12) 

where WAt[fc] is the logarithm of the stochastic variance. More specifically, we will choose ca* to 
be a gaussian white noise independent of w and of variance a 2 At. The choice for the process a; At 
introduced in Jl§f| , is dictated by the cascade picture. It corresponds to a gaussian stationary 
process whose covariance can be written 

Cov(w A t[fc],WAtW) = A 2 \np At [\k - l\] 
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where pAt is chosen in order to mimic the correlation structure observed in cascade models with 
an integral time T: 

f mrrjAi for |fe| < T/At — 1 



pAt[k] 



I otherwise 



Hereafter, we will refer to the process u>(t) as the "magnitude process". In order the variance of 
XAt(t) to converge when At — ► 0, one must choose the mean of the process to At such that [|l8| 

E(u! At [k}) = -Var(wAtW) = -A 2 ln(T/At), 
for which we find Var(XAt(i)) = o 2 t. Let us review the multifractal properties of MRW. 

4.2 Q spectrum: computation of the moments 

The qth-order moment of the increments of the MRW can be computed. Since, by construction, 
the increments of the model are stationary, the law of XAt(t + 1) — XAt{t) does not depend on t 
and is the same law as XAti})- In Ref. flgl , it is proven that the moments of X(l) = -Xa*— >-o+(0 
can be expressed as 

E(X(l) 2 n = ^^- [dm... [ du P Hp( Ui - Uj ) iX \ (13) 
P- Jo Jo t< j 

where p is defined by 

p(t) fr/|t| for|t|<T 
w [1 otherwise 

Using this expression in the above integral, a straightforward scaling argument leads to 

^ p-2p(p-l)A" 

T 

where we have denoted the prefactor 



M(2p,l) = K 2p [U , (141 



K 2p = TPa 2p (2p-l)U [ dux... [ du p T\ 

Jo Jo i<3 



Ui - u 3 1 



(15) 



Note that K 2p is nothing but the moment of order 2p of the random variable X(T) or equivalently 
of drXft). From the above expression, we thus obtain 

C 2p -p-2p(p-l)A 2 

and by analytical continuation, the corresponding full £ q spectrum is thus the parabola 

C - (q - q(q - 2)A 2 )/2 . (16) 

Let us remark that one can show that K q = +oo if ( q < (i.e., q > 2 + 1/A 2 ) and thus the 
pdf of 6iX(t) have fat tails fl8|| . In order to control the order of the first divergent moment 
(without changing A), one could simply choose for the eAt's a law with fat tails. Indeed, the 
prefactor a 2p (2p— 1)!! in Eq. ( fl5| ) comes directly from the fact that the eAt's have been chosen 
to be Gaussian. Using instead fat tail laws (e.g., t-student laws) would allow us to control the 
divergence of this prefactor. 

In Fig. H we have estimated the scaling behavior of the absolute moments M(q, I) for a 
discrete simulation of a MRW (Fig. ||(a)). In order to simulate the sampling of a time continuous 
MRW, we have generated a discretized MRW using At << 1 and then subsampled it at the 



9 



sample period 1. Using this procedure, we have generated a 2 17 long time-series using the 
parameters At = 1/16, T = 2 15 , a 2 = 1 and A 2 = 0.03. In Fig. |(b), we have plotted, in double 
logarithmic representation M(q, X) versus I for different values of q. In these representations, the 
linear behavior of each moment indicates that the scaling hypothesis is verified. The estimation 
of Q (made by estimating the slope of each of such curve) is reported in Fig |2|(c). As expected 
this spectrum is a parabola that is in very good agreement with expression ( jig ) . 

It is clear that the same power law scaling does not stand when I goes to +00. Since p(l) = 1 
for large I (as compared to T), we get 

CP 

Thus, there is a cross-over from the parabolic multifractal behavior at time scales I < T which 
is described by Eq. (|l^) to the Brownian-like behavior at larger time scales (I » T) 

Cg - 9/2 • 

In Eq. (^), we have shown that there exists a deep link between the self-similarity kernel 
and the £ q spectrum. This suggests that the probability distribution functions of our model 
satisfy Castaing's equation when going from large to small time scales with a gaussian kernel 
Gi x- Thus, as far as the increment pdf at different time scales are concerned, they will satisfy 
an evolution equation from "quasi-Gaussian" at very large scale (/ >> T) to fat tailed pdf's at 
small scales. This transformation of the pdf's is illustrated in Fig. [I|(a) where are plotted, in 
logarithmic scale, the standardized pdf's (the variance has been set to one) for different time 
scales in the range [1, AT) . The pdf's have been estimated for 500 realizations of size 2 17 of MRW 
with parameters A 2 = 0.03 and T = 2 13 . In solid line, we have superimposed the Castaing's 
transformation obtained from the coarse scale pdf (at scale T) using the appropriate normal 
self-similarity kernel. If Fig. |](b) we have reproduced similar analysis for the S&P500 future 
variations. Besides statistical convergence limitations, one can observe the same features as in 
Fig. |(a). 



4.3 Volatility and magnitude correlation functions 
4.3.1 Volatility correlation functions 

As recalled in the introduction, increments of financial time series are well known to be uncorre- 
cted (for time lags large enough) while their amplitude ( "local volatilities" ) possesses power-law 
correlations. Let us show that our model satisfies these two properties at all time scales smaller 
than the "integral time" T . By construction, the increment correlation function, 

((x At (t + t) - x At {t)){x At {t x + n) - x At (t x ))) 

(V \t\ — t\ > t), is zero in our model. Let us study the correlation function of the squared 
increments. Since the increments are stationary, we can choose arbitrarily t\ = 0. Thus we need 
to compute, in the limit At — > 0, the following correlation function, that corresponds to a lag I 
between increments of size r 

C(l, t) = ((X At (l + r) - X At (l)) 2 X At (r) 2 ). (17) 

From the results of Ref. |l^] and in the case < I < T, < r + I < T, we get, 

f l+T r 2 

C(l, t) = cr 4 / du dvp(u - v) 4X . 
Jl Jo 
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Figure 3: Multifractal Analysis of a MRW sample, (a) Plot of a sample time serie of length 
2 17 . The sampling size and the trend amplitude have been chosen arbitrarily to be compared to 
fig|l|(a). (b) Log-log plots of M(q,l) versus I for q = 1,2,3,4,5. The time scales I range from few 
minutes to one year, (c) ( q spectrum estimation (dots) and comparison to prediction as given by 
Eq. (|D (solid line). 



11 




Figure 4: Continuous deformation of increment pdf's across scales, (a) MRW Model. 
Standardized pdf's (in logarithmic scale) of 5iX(t) for 5 different time scales (from top to bottom), 
I = 16,128,2048,8192,32768. These pdf's have been estimated on 500 MRW realizations of 2 17 
sampled points with A 2 = 0.03, At = 1/16 and T = 8192. One can see the continuous deforma- 
tion and the appearance of fat tails when going from large to fine scales. In solid line, we have 
superimposed the deformation of the large scale pdf using Castaing's equation (||) with the normal 
self-similarity kernel. This provides an excellent fit of the data, (b) S&P 500 future. Standardized 
pdf's at scales (from top to bottom) I = 10, 40, 160 min, 1 day, 1 week and one month. As in Fig. 
(a) the scale is logarithmic and plots have been arbitrarily shifted along vertical axis for illustration 
purpose. Notwithstanding the small size of the statistical sample (as compared to (a)), one clearly 
sees the same phenomenon as for the MRW. The fact the Castaing's equation (||) allows one to 
describe the pdf's deformation across time scales of financial assets has originally been reported in 
Ref. H where similar plots for FX rates can be found. 
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A direct computation shows that 



r l+T a„. r v i-4a 2 



(1-4A 2 )(2-4A 2 ) 

and consequently 



Si T duJ^dv\u- 
((l + r) 2 - 4A2 + (I - r) 2 - 4A2 - 2; 2 - 4A2 ), 



C(l, r) = K(\l + r| 2 - 4A2 + \l - t| 2 - 4a2 - 2|^| 2 - 4a2 ) (18) 



where 

(J 4 T 4A 2 



K = 



(1 -4A 2 )(2-4A 2 ) 
Let us note that in the usual case < t << I, one gets 



I - - 4A 



C(Z,r)^ 4 r 2 ^-j (19) 
i.e., for fixed t, the volatility correlation function scales as 

c(/) - r 2v (20) 

with ^ = 2A 2 . From the estimates A 2 ~ 0.025 — 0.05 for financial assets (see section 5), one thus 
obtains v ~ 0.05 — 0.1, values very close to the ones observed empirically in many works. 



4.3.2 Power of returns and magnitude correlation functions 

Let us now show that magnitude correlation functions behave as expected, i.e, decrease very 
slowly as a logarithmic behavior. 

For that purpose, the previous computation of the correlation function can be extended to 
the power of returns \X At (l + T~) — X& t {T)\ p ■ Several empirical works have concerned the study of 
such "generalized volatilities" and people often noticed variations of amplitude of the correlation 
and of the power-law exponent v v when varying the order p J37| . In Ref. Jl8| , it is shown that 
the quantity, 

C p (1,t) = (\X At (l + T)-X At (l)\P\X At ( T )\P) , (21) 
behaves, when r is small enough, as 

— A 2 2 

C p {l,r)~Kl^) 2Cp (Ly P (22) 

where the constant K p has been defined previously. Using analytical continuation of the behavior 
of C p in the limit p = e — > 0, we can obtain, from previous expression, the behavior of the 
magnitude correlation function C^(1,t): 

C u (l, t) ~ 6" 2 (C £ (l, t) - M(e, r) 2 ) ~ -A 2 ln(A) . (23) 

The magnitude correlation function, for r small enough, has thus the same behavior as the 
correlation function of the underlying magnitude process u> A t- This result is checked in Fig. || 
where we have plotted the magnitude correlation function for r = 32At as a function of ln(Z). 
This correlation function has been estimated using a single realization of the process of 2 17 
sampled points, i.e, 16 integral scales. The linear behavior we obtain is exactly the same one 
as predicted from Eq. ( p3| ) and Fig. [| Measures of the slope and the intercept of such straight 
line provide a good estimate of respectively A 2 and T. 
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Figure 5: Magnitude correlation function of the model. The correlation function has been 
estimated on 16 integral scales (see text). In continuous line we have superimposed the correlation 
function of the magnitude process WAt that is involved in the stochastic volatility. 

4.4 Extension to a multivariate process 

In order to account for the fluctuations of financial portfolios and to consider management 
applications of our approach, it is important to build a multivariate, i.e. a vector valued, 
version of the previous multifractal random walk. Since only gaussian random variables are 
involved in the construction of section 4, this generalization can be done by considering two 
uncorrelated gaussian random vectors eAt(t) and LUAt(t) whose covariance matrices are denoted 
respectively £ and A. Hereafter, we will refer to these matrices as respectively the "Markowitz 
matrix" S and the "multifractal matrix" A. One can then define the multivariate multifractal 
random walk (MMRW) X(t) as: 



with Cov(ei(t),ej(t + t)) = <5(r)Sy and Cov(wj(i), Uj(t + r)) = A+j ln(Tij/\At + |r|) (note 
that the previously defined coefficients a 2 and A 2 for an asset i correspond respectively to the 
diagonal elements and An). Let us briefly review some of the properties of this model, 
postponing its detailed analysis to a forthcoming publication |l7[ . A quantity that will be of 
central interest is the fc-point joint moment of order q\, q 2l ■■■,qk that can be defined as: 



Xi(t + At) - Xi(t) = 6i(t)e^ , 



(24) 



M ll _ lk (q 1: ...,q k )^E(\X i M qi ---\X lk (l)n . 



(25) 



When k = 2, by denoting i\ = i, %i = j, q\ = p and qi = q, let us define the joint scaling 
exponent spectrum as: 




dj (pj q) = d (p) + (?) - A ijpq , 



(27) 



where Ci(#) is the £ 9 spectrum for the component Xi(t). Thus, for uncorrelated oVs, one has 
CijiPtQ) = dip) + Cj(q) while for the extreme case u>i = u>j, the exponent becomes (ij{p,q) = 
Ci(p + q) — Cj(p + q)- The computation of the scaling exponent is trickier for general Markowitz 
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and multifractal matrices. Under some mild conditions that are necessary for the existence of 
a non trivial limit At — > + , one can show that the previous scaling law remains valid even for 
non diagonal matrix X ]rj| . 

In order to define a simple way to get an estimate of the multifractal covariance coefficient 
Ay, let us define the moment ratio: 

From Eq. (^) , the value of Kij (q) is simply 

K i:i (q) = ~A ljq 2 . (29) 

Thus, the non-diagonal clement in the multifractal matrix A corrresponds to the nonlinear 
behavior of the exponent spectrum n{g) of the moment ratio R. Along the same line as for the 
computation of the magnitude auto-correlation in previous section, one can get the correlation 
function of magnitudes LUi(t, I) and u>j(t, I) from the limit q — * of Rij(q, I): 

Cov(^(i,0,^(U)) ~ -AijMO+C, (30) 

where C is a constant related to T^- [[l7| . Thus the scale behavior of the magnitude covariance 
provides an estimate of the multifractal correlation coefficient Ay . This is the generalization of 
the classical result in multifractal analysis that relates the intcrmittency coefficient A 2 = An to 
the scale behavior of the variance of the magnitude. 

Let us remark that the covariance of the variations of the assets i and j can be obtained by 
a direct calculation: 

Cov(Xi(l), X 3 (l)) = Sj . e ^«+A 3 - 3 -2A y )j _ (31) 

This covariance between Xi(t) and Xj(t) thus depends not only on S, the "Markowitz" co- 
variance matrix, but also on the multifractal matrix A. This expression, allows us to get an 
estimate of the value of once the values of A are known. 

Finally, let us mention that the idea of "multivariate multifractality" has been recently intro- 
duced in Ref. where the authors propose a phenomenological generalization of Castaing's 
equation to the multivariate setting. Evidences that financial assets are characterized by non 
trivial multifractal matrices are also provided. We are currently working to obtain further em- 
pirical evidences towards such conclusions. Moreover, a precise link between the present model 
and the extended Castaing's approach of Ref. is under progress. 



5 Parameter estimation for real financial data 

We have seen that the MRW is characterized mainly by 3 parameters: cr 2 , the white noise 
variance, T the integral scale and A 2 the magnitude variance. We have shown that this model 
is able to reproduce all the main features of the future S&P 500 time series. Natural estimators 
of those parameters can be defined from the results of previous section. The parameter A 2 
can be obtained from the shape of the £ q spectrum that is itself estimated using the scaling 
behavior of the absolute moments M(q,l). This parameter can also be estimated thanks to the 
magnitude correlation function 0^(1, r) that behaves as —A 2 ln(Z/T). From the intercept of such 
correlation function as a function of m(Z), we can define an estimator of the integral scale T. 
Finally, the parameter a 2 can be obtained using the classical relationship Vax(SiX / \t(t)) = <r 2 l. 
In this section, we report estimates of the multifractal parameters A 2 and T for some financial 
time series. We do not have the ambition to provide fine estimates of those parameters. Our 
aim is rather to get an idea of realistical values of the parameters of the model for real assets. A 
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Series 


Size 


X 2 


T 




7 1 n 4 


U.U AO 


3 yGcirs 


Future TY/TTSD 


7 1 n 4 


09 


ri TDAnt nQ 


Future Nikkei 


7.10 4 


0.02 


6 months 


Future FTSE100 


7.10 4 


0.02 


1 year 


S&P500 index 


6.10 a 


0.024 


3 years 


French index 


6.10 a 


0.029 


2 years 


Italian index 


6.10 3 


0.029 


2 years 


Canadian index 


6.10 3 


0.024 


3 years 


German index 


6.10 a 


0.027 


3 years 


UK index 


6.10 3 


0.026 


6 years 


hong-kong index 


6.10 a 


0.05 


3 years 



Table 1: Multifractal paramater estimates for various assets 



precise discussion of the properties of various estimators from a statistical point of view is out 
of the scope of this paper and will be addressed in a forthcoming publication. Note that similar 
empirical study has already been performed in Ref. [ fl9| . We have studied some high frequency 
future time series that are sampled at a 10 min rate over the 7 years period from 1991 to 1997. 
We have also processed a set of daily index values for 8 different countries over the period from 
1973 to 1997. The results are reported in table [l]. 

We remark that the values of the multifractal parameter A 2 are all very close to 2.5 10 -2 
(excepted for the hong-kong index) . The integral time T values are centered around 3 years but 
with a large spread. Let us notice that we get an estimate of hi(T') and thus the error on the 
estimate of T can be very large. We do not report here the values of the errors and confidence 
intervals for the proposed estimators that will be studied elsewhere. 

6 Discussion about other approaches and findings 

In this section, we make some comments about related studies that concern multifractals and 
finance. 

6.1 Turbulence and finance 

The analogy between turbulence and finance has been originally proposed by Ghashghaie et 
al. These authors proposed to describe the pdf's of FX price changes at different time 

scales in the same way physicists describe the pdf's of velocity variations at different space 
separations in fully developed turbulence. This approach naturally leads to the notions of 
cascading process, Castaing's formula and multifractality as described in section 2. This work 
suggests that the key mechanism at the origin of these observations, is an information cascade 
according to which short-term traders are influenced by long-term traders. This cascade is the 
analog of the Richardson's kinetic energy cascade in turbulence where small eddies result from 
the breakdown of larger ones and so on |H] . If the observations reported in Ref. Jl(J strongly 
support this point of view, its quantitative understanding in terms of "microscopic" mechanisms 
remains an open question. In this section we would like to comment about some criticisms that 
have been raised about the analogy between turbulence and finance. The first one concerns 
the power spectrum behavior in both situations jn], |lj, |). In turbulence, Kolmogorov theory 
predicts a fc -5 / 3 power spectrum that is confirmed in experimental situations. In finance, since 
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price fluctuations are almost uncorrelated, they are characterized by a k~ 2 spectrum. For a 
general multifractal process, the exponent (3 of the power spectrum behavior can be shown to 
be related [^7], |2^| to the value of £2: /3 = 1 + £2- Thus, from the cascading process point of view, 
nothing prevents the exponent /3 from being equal to the exponent of the Brownian motion, i.e., 
= 2. In other words, as examplified by the MRW, a cascading process can have uncorrelated 
increments. We could also remark, that in turbulence (3 = 5/3 has a dimensional origin, i.e., it 
is the exponent of the spatial spectrum of velocity fluctuations within an Eulerian description. 
If one adopts a Lagrangian description and one is interested by temporal fluctuations of a fluid 
particle velocity, then the dimensional value of the power spectrum exponent is (3 = 2. Thus 
the value of this exponent is not a pertinent argument to reject the analogy with turbulence. 
Another difference that has been raised in [Q] concerns the behavior of the probability of return 
to origin P;(0) that has been shown to possess a scaling regime in finance while its behavior 
is more complex for a turbulent velocity field. First of all, let us point out that whatever the 
quantity studied (probability of return or absolute moments), it is well known that there is no 
observed well-defined scaling regime in turbulence: the classical "log-log" plots always display 
some curvature across scales. This curvature is Reynolds number dependent and several studies 
suggest that it vanishes, i.e. the field is scale-invariant, only in the limit of infinite Reynolds 
number |2l], |3l], [3S(] . However, within the cascade paradigm and using Castaing's equation, the 
scaling behavior of the probability of return to origin is easy to show. Indeed, by setting SX = 
in ([5]), one obtains, from the definition of £ 9 and the self-similarity kernel: 



The exponent for the probability of return to origin is thus simply C—i- For the log-normal 
stochastic volatility model introduced in section 4, we thus get 



To conclude, neither the power spectrum exponent, nor the scaling behavior of the probability 
of return to origin can be used as argument against the existence of a cascading process at the 
origin of the fluctuations of financial time series. 

6.2 Subordinated processes. Multifractal time 

Subordinated processes are Markov processes in a time variable ^i(t) that is itself an (increasing) 
random process [ |30| . Such processes have been introduced in finance by Mandelbrot and Taylor 
p6{ to account for the existence of Levy stable laws as the result of a Brownian motion in 
some stochastic time. Today, the idea of modelling financial return fluctuations as a Brownian 
motion in a "fractal time, "trading time" or "financial time" can be found in many approaches. 
In Refs. |r| [l4|, the multifractal nature of these fluctuations has been modelled by a (fractional) 
Brownian motion subordinated with a multifractal stochastic measure. In this section, without 
any concern for rigor, we would like to make a link between our stochastic volatility approach and 
the multifractal time approach of Mandelbrot and co-authors. Let us first remark that if we drop 
the noise e in Eq. ([l2] ) and keep only the stochastic volatility a(t), we can construct a stochastic 
measure fi(dt) that satisfies [i(dt) = e^&'dt. Using exactly the same kind of computation as in 
section 4, one can show that this measure is stationary and its multifractal spectrum r(g) is 




(32) 



(33) 



r(q) = C2g - 1 , 



(34) 



as usually defined by 



M[cm])«>~^) +1 



(35) 
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Let us note that the existence and the construction of such a measure that is stationary and 
possesses a continuous scale invariance, was at the heart of the construction in Refs. |ll| [l4| 
and was still an open problem. According to these studies, one can thus construct a multifractal 
process by simply considering the subordinated process S(t) — £?(/i([0, i])) where B(t) is the 
standard Brownian motion. The £ g spectrum of such process would be exactly the same as 
the stochastic volatility process defined in section 4. This is not so surprising since, formally, 
a differential form for the process S(t) would be dS = ^dB{^{t)). If one assumes that a 
white noise that is subordinated remains a white noise, one thus obtains dS — e u ^'dB{t) that 
is the equation that defines the MRW of section 4. The questions of well-definiteness of this 
construction, its statistical properties and the precise mathematical justification of such results, 
will be addressed in a forthcoming work. 



6.3 Some remarks about Bouchaud, Potters and Meyer's model 

Besides multifractal and cascade pictures, our present approach has been inspired by a recent 
paper by Bouchaud, Potters and Meyer |3q] . These authors have proposed a model that is very 
similar to ours: the stochastic volatility <x(i) instead of being log-normal (e w W) is a normal (u>(i)) 
random process with long-range (power-law) correlations. By a simple analytical computation, 
they have shown that the q-order cumulants of such a process satisfy a simple scaling behavior 
but the moments display apparent multiscaling caused by a "competition" between the different 
cumulant behavior on a finite scale range. They thus conclude that a distinction between 
multifractality and such "apparent multifractality" is a difficult task for finite size time series. As 
far as multifractal analysis and modelling of financial time series are concerned, this work is very 
interesting and the previous assertion is undoubtedly difficult to infirm. However, let us remark 
that in order to illustrate their purpose, Bouchaud et al. choose a "stochastic volatility" <r(t) — 
gw(t) i ns tead of their "monofractal" model a(t) = \u)(t)\. The reason invoked by the authors is 
that the log- normal is "a more realistic time series as compared with real data... without changing 
the feature of the above model, i.e. the very slow decay of the volatility correlations" . They thus 
claim that the scaling features of both models are the same and thus that the multifractality 
observed for the simulations of the "log-normal" volatility model is only apparent as predicted 
by their theory for the "normal" volatility model. The results reported in section 4 can be used 
to show that this interpretation is not correct. Let us indeed reconsider both results of Ref. 

and section 4. According to the "normal" volatility model, the moment of order q = 2p is 
written in terms of cumulants and behaves as ps| : 

M(2p,l = NAt) = A 2p , N^-^ p + ... + A 2 ,.. }2 NP (36) 

where the constants Aq u q h depend only q^ A 2 the variance of u> and v the exponent for the 
correlation function of u>: C W (Z) ~ l~ v . According to this equation, if N is small enough, 
M(2p, I) ~ lO--")P while, for N very large, M(2p, I) ~ P. The transition scale N*(q = 2p) above 
which the scaling exponent is £ g = q/2 can be estimated if we define it as the scale where the 
contribution to the moment of order q of the cumulant of order 4 and 2 are equal. Using the 
expression in Ref. ps| for second and fourth cumulants, C 2 and C4, we can show that at scale 
N*(q = 2p), we have (2p ~ l)!!Cf ~ p(p - l)(2p - l)!!Cf ~ 2 C* 4 /6. From the value of C* 4 , we 
obtain (q > I): 

N*{q = 2p) = p(p - l)i/ a 2 2 ("- 1 ) J2 ™ 2( "~ 1} (37) 

\ m=l / 

This function only depends on q and v and is increasing as q — » +00. Thus, the larger the q value, 
the wider the range of scales on which apparent multifractality exists. However, a numerical 
computation of the values of N* for v = 0.2 shows that the value N* ~ 100 is reached only 
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for the moment of order p = 15. The greatest moment value attained in practical situations is 
q ~ 6, for which N* < 1 ! That means that for all moments less that 10, the model of Bouchaud 
Potters and Meyer predicts the trivial spectrum £ q = q/2 without any cross-over phenomenon. 
For their numerical simulations, they have used a log-normal model. However, within the log- 
normal ansatz, the conclusions of Ref. p8j are questionnable since, when v is small enough, this 
model is very close to the model introduced in section 4. Let us indeed consider as in j38) that 

C«(0 ~ A 2 (r(^) cos(^r)) 2 l~ v with v very smallQ By expanding this expression, we obtain 

A 2 

C u (l) = A 2 ln(Z) + 0(fln(Z)) . (38) 

v 

If we set T — e v , then for 1 < I « T, this equation becomes 

C u (!) = A 2 ln(T/0 (39) 

that is the same correlation function as introduced for the multifractal model in section 4. Let 
us notice that for v = 0.1 we have T w 2.10 4 , T w 3.10 5 for v = 0.08 and T w 5.10 8 for v = 0.05 
!. In this model T is increasing very fast as v goes to zero. We can thus conclude, that the model 
numerically studied in Ref. fl38| can be seen as multifractal from one point of view: whatever 
that scaling range [1,T], there exists v small enough {y k, 1/ ln(T)) such that the model displays 
multiscaling with log-normal ( q spectrum in this scale range. According to these remarks, we 
thus think that the multifractal picture is more realistic to describe multiscaling in financial 
time series. 



7 Summary and prospects 

In this paper we have reviewed what are the main features of multifractal processes. We have 
shown that the Multifractal Random Walk is a very attractive alternative to classical cascade 
processes in the sense that it is stationary, continuously scale-invariant and formulated using a 
simple stochastic evolution equation. As a model for financial engineering, MRW are interesting 
for many reasons. First, as illustrated in details for the S&P 500 intraday time series, this model 
is able to reproduce the main empirical properties observed for financial time series. Moreover, 
as Brownian motion and other stable walks, it is a "scale-free" model in the sense that it does not 
have to fit a particular time-scale since it is scale-invariant. This kind of stability with respect 
to time "aggregation" is a serious advantage as compared to classical ARCH-like models which 
parameters strongly depend on the time-scale one is interested in. Moreover, as discussed in 
section 4.4, a simple multivariate formulation of MRW can be proposed. To our knowledge, it is 
the first example of an extension of the notions of multifractality to a vector field. The empirical 
results reported in Ref. |l^] suggest that MMRW can be pertinent for portfolio theory. We are 
currently working on applications of MRW to classical problems of finance like management 
problems and option pricing theory. 

From a theoretical point of view, MRW can be seen as the simplest model that contains the 
main ingredients for multifractality. In that respect, it can be very helpful to elucidate, in many 
fields where multiscaling is observed, what are the generic mechanisms that are involved leading 
to "non-trivial" self-similarity properties. Various "microscopic" models, as proposed in finance 
or other fields, could be considered within this perspective. It could also be interesting to recast 
our approach within a field theoretical formulation involving some renormalization procedure. 

7 Notice that there is no reason to consider the same value of v for the normal and log-normal models. The results 
of this paper suggest that, in finance, the value for the correlation exponent in the log-normal model is very close to 
zero and significantly smaller than 0.2 
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From a mathematical point of view, this problem is deeply linked to the existence of a limit 
stochastic process when the sampling time At goes to zero. The convergence of the moments is 
not sufficient to prove this non trivial assertion. Such a limit could be very useful to develop a 
new stochastic calculus within which, for example, one could formulate the model of multifractal 
time of Mandelbrot and co-authors very naturally (see section 6.2). Finally, in a forthcoming 
work, we will discuss the generalization of such approach to other laws than the (log-)normal. 
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